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Abstract 

A higher spin generalization of the XYZ spin chain defined by the 
fusion procedure is considered. Energy and momentum of the ground 
state and low-lying excited states of the model are computed by means of 
an isomorphism between a space of theta functions on which the Sklyanin 
algebra acts and a space obtained by the fusion procedure. 
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1 Introduction. 



In this letter we calculate energy and momentum of low-lying excited states of 
the higher spin generalization of the XYZ spin chain model by means of the 
algebraic Bethe Ansatz method. This work is a continuation of the author's 
previous works 0] , || , where a generalization of the eight- vertex model || was 
studied and the two particle S'-matrix of the corresponding spin chain model 
was computed. There we assigned a (21 + l)-dimensional space to the vertical 
edges of the lattice and a two-dimensional space to the horizontal edges. In 
the present paper we make use of the fusion procedure (see || , || , || , ) to 
assign a higher dimensional space to the horizontal edges and then consider the 
corresponding one-dimensional quantum spin chain model. 

We shall give in an explicit isomorphism of a representation space of the 
Sklyanin algebra which is defined as a space of theta functions with the space 
of symmetric tensors. (This is a special case of more general isomorphisms by 
Hasegawa ||.) Through this identification we can identify the L operator of the 
model in || | with a special case of R matrices in Q , and use the results in to 
compute the general models. In particular this identification is indispensable to 
interpret a special value of the logarithm of a transfer matrix and its logarithmic 
derivative as a momentum and an energy operator. See, e.g., 

As is conjectured in B from the corresponding results for the higher spin 
XXX and XXZ models (l0| , , jl2| , the energy and the momentum are inde- 
pendent of the spin of the local quantum space and expressed as a sum of two 
terms each of which depends on a rapidity of a hole in the string configuration. 
Thus we can justify the interpretation of those states as "two-particle states" . 

We make use of the results in Q] , but adopt simpler normalizations in Jl3[ . In 
particular, we normalize the R matrices by the unitarity condition, and therefore 
they are meromorphic functions of the spectral parameter, while || and fl^| use 
holomorphic R matrices. 



2 Fusion procedure. 

In this section, we briefly review the fusion procedure for the elliptic R matrices. 

Let V 1 = Sym(Vi <8> • • • ® V 2t ) for I S |Z, where V- C 2 (i = 1, . . . , 21) 
and Sym is the symmetrizer. The elliptic R matrices R 1,1 (u) have the following 
properties. 

(i) R 1 ' 1 (u) is a linear endomorphism of V 1 (g> V 1 meromorphically depending 
on a complex parameter it. 

(ii) ( Yang-Baxter equation) As an endomorphism of V ® V 1 ®V l , 

i?12 (ui - u 2 )R l i3 (ui - u 3 )R l 2 ^ (u 2 - u 3 ) = 

= i?23 ( (u 2 - u 3 )R l { l 3 (ui - u 3 )R l { 2 (ui - u 2 ). (2.1) 
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(iii) ( Unitarity) As an endomorphism of V <8> V , 

i?y (u - w)i?2i'(w ~u) = Id vl ® vl > . (2.2) 

(iv) When u = 0, R l ' l (u) is a permutation operator: for all v,w E V , 

#*''(())(« ® to) = to<g>u. (2.3) 

They are constructed in || , Q , by the fusion procedure [Q from Baxter 's 
i? matrix i? 1 / 2,1 / 2 (u) = i?(w; r) defined by 



o=0 



2^(r ? ;r)^ 1 (u + 2r ? ;r)' (24) 



where go = (H), gi = (10), 52 = (00), g$ = (01). The explicit definition of R 
is: 

R^'Au) :=Sym w J R VitV _(u+(2l'-l)rj)--- (2.5) 
■ • • R Vt .v 7 (u + (2j - 21' - l)r;) • ■ ■ R Vi ,v r (u + (-21' + 1)77), 

(2.6) 

R l / l vl , (u) := Sy mi ... 2i tf££ v (u + (21 - ■ ■ ■ (2.7) 

• • • < /2 ^, (u + (2j -2i- 1)7,) . . . i?^;, ( M + (-21 + 1),?). 

(2.8) 

Here Vi = Vj = C 2 , suffixes of R designate the spaces on which the R matrix 
acts and Symj m is the symmetrizer on the space V\ (8> ■ • • <8> V m etc. 

There is another expression of R 1 ^ 2 ' 1 (u) in terms of the representation of the 
Sklyanin algebra which we used in ^). The Sklyanin algebra [Q U T . r] (sl(2)) is 
generated by four generators S°, S 1 , S 2 and S 3 satisfying the relations coming 
from the Yang-Baxter type relation of the L operator, L(u), defined by 



L(.)=5>>K^, %™ (2.9) 



,(1 

The Sklyanin algebra has a representation on a space of theta functions 

: f/ T ,,(s/(2)) -> Endc(eft + ), (2.10) 
f(y + l) = f(-y) = f(y), \ 



6qo = s /(y): holomorphic on 



/(2/ + T)=e- 4M (2y+r) /(2/) /- (2 n) 



It is easy to see that dim Oq[) + = 2/ + 1 . The generators <S a act on this space 
as difference operators. See |lj| and Appendix A. of 0], |l3|]. We fix a basis of 
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ejj+ (9 00 (2y,2T) 



10 (2y; 2r), 6> 00 (2y; 2r) + lo (2y; 2r)) and identify 6 2 + with 
C 2 through this basis. Then we fix an isomorphism of the space of symmetric 
tensors and the spin I representation space of the Sklyanin algebra as follows: 

V 1 = Sym(Vi ® • •• <g> V21) 9 Sym(/i(yi) ® • •• ® fiiiyii)) >-> 



h{y)---hi{v) cv l = e 



4H 



(2.12) 



where = C 2 (i = 1, ... ,21) identified with OqJ , fi(yi) S V; and Sym is the 
symmetrizer. Under this identification the L operator ( |2.9| ) is proportional to 
the i? matrix defined by (2.6): 

6 n (2r)) 



R^iu) 



Id C 2 ®pW(L(u + rj)). 



(2.13) 



? 11 (u+(2/ + l)r ? ) 

This can be verified by comparing the action of both hand sides on the inter- 
twining vectors. (See Lemma 2.1.3 and Theorem 2.3.3 of |(| (Adv. Stud. Pure 
Math. 16) for the left hand side and (1.18-21) of (ll| for the right hand side. 
See also Q for general cases.) 

As is the case with the trigonometric and rational R matrix, there is a 
recurrence relation with respect to the auxiliary spin |[(|, p2[] . Let R 1,1 (u) and 
R 1 / 2 ' 1 (u) are the R matrices on the space Sym(^(8)- ■ • V\)®V l and on Vq®V 1 
respectively. Here Vi = C 2 (i = 1, . . . , 21) and V 1 is a space of symmetric tensors 
defined above. Then as an operator on Sym(V2z ® • • • Vi) ® Vq ® V , 



i? M '(u + ?7)i? 1/2 ^'(i 



- 2/77) = 



(«) 



q-detE 1 / 2 - r (u-(2Z-l)?7)x 
xR l - 1 ^ l '(u + 2r 1 ) 



(2.14) 



where q-det i? is the quantum determinant p7[ defined by 
q-deti? 1 / 2 ' i '(u)=troiP ii? 1/2/ (u + ?7) J R 1 / 2 '''(u-r ? ) = ^ 



i(u-2l'r)) 



viu + 21'^) 



Id 



(2.15) 



(P ~^ is a projection to the antisymmetric tensor in Vi (8> Vo). The block structure 
of the right hand side of (2.14) comes from the decomposition of the tensor 
product of the auxiliary spaces by the Young symmetrizcrs as follows: 



1 


2 




21 


<8> 










1 


2 




21 


© 


1 


2 




21 
















(2.16) 



Here we denote the image of the Young symmetrizer by the corresponding Young 
tableau. Since the second Young symmetrizer in the right hand side of (2.16) 
gives an isomorphism from the space Sym(Va; ® • • • <8 V2) (8 Ant(Vi <g> Vq) = 
Sym(V2/ ® • • • <S> V2) to its image, we identify these spaces in the right lower 
corner of the right hand side of (2.14). The proof of this equation is based on 
Lemma 5 of pi and done by the same argument as that in section 1 of p6[. 
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3 Higher spin XYZ model. 

In this section we define a higher spin generalization of the XYZ model and 
apply the algebraic Bethe Ansatz. 
The state space of our model is 

H = V±®V±®---®Vj r , (3.1) 

where Vj = V . We define the model by specifying the transfer matrix, namely, 
generating function of quantum integrals of motion, as follows: 

T l '' l (u) := tr vl , Ryj, («) ■ ■ • R 1 ^, («) G End c (W), (3.2) 

' N ' 1 

where R v i> v j. acts non-trivially only on the component V 1 V~ of V 1 <g) H. 
The most important property of the transfer matrix is the commutativity 

[T l '> l (u),T l "> l (v)} = 0, (3.3) 

which is a cons equ ence of the Yang-Baxter equation (2.1). 

Thanks to (2.3), we can define a momentum operator p and a Hamiltonian 
H of the spin chain by: 



p = * logT M (0), H = const. ^-\ogT l ' l (u) 
1 ciu 



(3.4) 

U=0 



where const, is a constant which we do not fix here. (See ||, |Q, jnj. |TT|. 

Due to the recurrence relation of the R matrix (2.14), the transfer matrix 
satisfies the following recurrence relation: 

= T^u + ,)T^\u - 21'rj) - ^ ( ;+^;;+ 1 7^ T''-^( M + 2,). 

0ii(it+(-2<' + l + 2Z)7j) ^ 3 5 ^ 

Hence the diagonalization problem of T l ' l (u) reduces to that of T 1 / 2 ' 1 ^). 

Hereafter we assume that the elliptic modulus r is a pure imaginary number 
r = i/t, t > and the anisotropy parameter r\ = r' jr is a rational number. We 
also assume that M = Nl is an integer. 

It was shown in jj| that there exist vectors & v (wi, . . . ,wm) depending on 
an integer v and complex parameters (wi, ... , wm) such that 

T 1 / 2 ' l (u) i & l ,(wi, . . . ,w M ) = t 1/2 < l (u;v,w 1 ,... ,w M )^u{w u ... ,%), 

(3.6) 



4 



provided that (y;wi, . . . , wm) satisfy the Bethe equations: 



hijwj +21ti;t) \ _ 4nivvTJ 0u (wj-w k +2r); t) 



N M 

0\\{w 3 -2lr);r) J On(wj - Wk - 2t\\t) ' 

for all j = 1, ... , M. The eigenvalue t 1 / 2 ''^) is equal to 

«,„) - + *(«)%T^ (3-8) 

where 

Q(u) := e g en(u - w 3 + V ), h(u) := { —-^^ j 

Inductively using the recurrence relation (3.5), we can prove that the eigenvalue 
of T l '' l (u) is 

t • l (u; i/,wi,... , wm) = 

= V j'.i, v Q(" + (2*' + 1)t?)Q(u - (21' + 1)77) 

^ j [ ) Q(u+(2l' + l-2j) V )Q(u+(2l'-l-2j)r 1 y l ' j 

where a l ^ l {u) = 1, ^'''(u) = Ul=i K u + ( 2// + 1 - 2k )v)- 

4 Thermodynamic limit. 

In this section, making use of the results of || , we compute several thermody- 
namic quantities of our spin chains. 



First let us recall several facts on the solutions of the Bethe equations (3.7) 
that we found i n |2j| . These solutions satisfy the string hypothesis which goes 
back to Bethe |19| . For later convenience we rescale the parameters as fol- 
lows: x = itu, Xj = itwj. In this notation the string hypothesis says that for 



sufficiently large N solutions of (3/7) cluster into groups known as A-strings, 



A = 1,2,..., with parity ± and centre xf'^'- 

= 4* + ^ + 0(e" W ), ° = Z 4^> ■ • • , ^ 

1 1 1 (4.1) 

where Imcc^' + = and Imiy = t/2. We denote the number of A-strings with 
parity ± by $(A,±). 

We consider the following string configurations which are consistent with the 
constraints found in p0[ when r\ = r'/r, r, r' are integers mutually coprime, r 
is even, r' is odd, and 2(2? + 1)t] < 1, as assumed in M. 
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• Ground state: u = 0, |(2Z, +) = iV/2, (t(A,±) = (t(2/, -) = for A f 21 
and centres of 2/-strings distribute symmetrically around 0. 

• Excited state I: (t(2i, +) = N/2 - 2, §(2/ - 1, +) = 1, |j(2Z + 1, +) = 1. 
There are two holes in the distribution of 2Z-strings which are denoted by 
x± and X2 and regarded as continuous parameters of the configuration. The 
Bethe equations determine the coordinates of the centres of the (21 ± lu- 
string which are denoted by x±. There are two possibilities: (x_,x + ) = 
((x 1 +x 2 )/2,(x 1 +x 2 )/2) or (x_,s+) = ((x 1 + x 2 )/2,( Xl + x 2 + l)/2). 

• Excited state II: tf(2Z, +) = N/2 - 1, Jj(2Z - 1, +) = 1, (t(l, -) = 1- 

There are again two holes in the distribution of 2^-strings which are denoted 
by x\ and x 2 . The Bethe equations determine the coordinates of the centres 
of the (21 — l)-string and the 1-string with parity — which are denoted by 
X- and xo respectively: (x~,xq) = ((x\ +x 2 )/2, (x\ + x 2 )/2) or (x_,a;o) = 
((xi+x 2 )/2,(xi+x 2 + l)/2). 

Using the results of [0 , we obtain the following asymptotics of the eigenvalue 
of the transfer matrix tr l (u) for large N. The largest eigenvalue of t 1 ' 1 which 
corresponds to the ground state is: 

^-\ogt l ' l (u; ground state ) ~ 
i 

»t / , „ , ., ,-. v~"^ sinh47rnZ7?f sinh7rni(l — Alif) . „ \ 

~ iV 7rZ + 2ttZx(1 - 4/77) + 2 > — — — l - sin 2nnx . 

y n sinh nnt sinh innrft X4 2) 

The excited states I and II have the same eigenvalue of the transfer matrix 
t l < l (u). Namely, 

— log t l,l (u; excited state I or II) — r \ogt l ' l (u; ground state) 
i i 

= logr(a; - x\) + logr(x - x 2 ), (4.3) 

where 

tt ^ sin27rna: 

logrfa;) := -- -ixx — > — . (4.4) 

' 2 ^ ncosh2?r™t v ; 

n—l 



(see |2l[ and Pj for details of computations). By the definition (3.4), the mo- 
mentum and the energy of these excited states are expressed as p = p\ + p 2 , 
H = H 1 + H 2 where 

^ l 2 Xl n cosh 2irnnb ' 

n—l 

I cos2irnxi \ 

Hi := const, -tt - 2tt > . (4.6) 

I ^ cosh27TO77i J v ; 



6 



Thus they are regarded as two particle spin wave modes each particle of which 
has rapidity Xi, momentum pi and energy Hi. In particular (4.5) and (4.6) show 
that the dispersion relation of these particles do not depend on I. 
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